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ABSTRACT 

1  n  22 

From  the  reaction  F  ^(a,  ny)Na  ,  the  low  excited 
22 

states  of  Na  have  been  studied  using  the  method  of  triple 
angular  correlations. 

At  a  bombarding  energy  of  5*38  MeV,  at  which  the 

yield  of  the  radiation  reaches  a  maximum,  neutron  gamma- 

ray  correlation  and  gamma-ray  distribution  experiments  were 

22 

performed  for  the  587  kev  and  the  892  kev  levels  of  Na 
The  beam-current  was  limited  to  a  maximum  value  of  0.2 
microamperes.  Because  of  the  low-coincidence  counting  rate, 
the  investigation  of  the  1.532  MeV  level  has  not  been 
performed . 

The  previous  probable  assignment,  1+,  for  the 
892  kev  level  does  not  fit  the  experimental  results  which  in¬ 
dicate  rather  a  2+  assignment.  The  two  unquantized  parameters 

of  the  angular  correlation  formula  have  been  measured.  The 

2  p  ( 1 ) 

population  parameters  ratio,  X  =  p-Kyy,  is  0.20  ±  0.20 
and  the  quadrupole  to  dipole  mixing  ratio,  6  is  equal  to 

0.72  ±  0.31. 
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CHAPTER  I 


1.1  Introduction 

The  knowledge  of  spectroscopic  parameters  like  spin 
and  parity  of  excited  states  or  the  multipole  mixing  ratio 
of  decay  radiation  is  very  important  for  a  better  understanding 
of  nuclear  matter  and  more  particularly  of  nuclear  forces 
which  are  revealed  indirectly  through  these  parameters.  Even- 
even  nuclei  have  been  investigated  more  and  more  in  recent 
years,  nuclear  models  describing  these  nuclei  have  been 
consequently  refined  and  the  matching  between  theoretical 
predictions  and  experimental  results  is  improving.  On  the 
other  hand,  nuclei  whose  configuration  is  far  from  closed 
shells  such  as  odd-odd  nuclei  are  more  difficult  to  compare 
with  models  proposed  by  the  theoreticians.  In  a  sense,  this 
should  excite  the  interest  of  experimentalists  because  some 
kind  of  "unified"  model,  allowing  an  understanding  of  all 
nuclear  forces,  may  be  expected  when  the  amount  of  information 
garnered  from  all  nuclei  becomes  sufficient. 

The  first  aim  of  this  experiment  was  to  determine 
the  spectroscopic  parameters  of  the  levels  892  kev  and  1.532 
Mev  of  Na  from  the  reaction  F  ^(a,  ny)Na  by  the  use  of 
triple  angular  correlations  as  described  by  Litherland  and 
Ferguson  (Li  6l).  For  reasons  explained  in  detail  later  (sec¬ 
tion  2.2),  the  investigation  of  the  1.532  Mev  state  had  to  be 


given  up. 


1 


2 

The  originality  of  the  method  proposed  by  Lither- 
land  and  Ferguson  lies  in  the  fact  that  if  a  nuclear  state 
is  formed  by  absorption  of  unpolarized  particles  in  the 
direction  of  the  quantization  axis  followed  by  the  emission 
of  a  second  particle  which  is  detected  along  this  axis,  then 
the  magnetic  quantum  numbers  of  the  substates  which  can  be 
populated  do  not  exceed  the  sum  of  the  spins  of  the  target 
nucleus  and  the  incident  and  emergent  particles. 

In  cases  where  the  sum  of  spins  does  not  exceed 

1 8  of 

for  instance  in  a  reaction  such  as  <j  (a,  n7)Ne^  ,  the 

angular  correlation  of  gamma  rays  emitted  by  deexcitation 
21  * 

of  Ne  is  independent  of  the  details  of  formation  of  the 
intermediate  state  and,  for  a  given  spin  assignment,  this 
angular  correlation  is  a  function  of  one  parameter  only, 
namely  the  multipole  mixing  ratio  6.  Since  angular  correlation 
experiments  generally  yield  the  two  values  aQ/aQ  and  a^/aQ, 
this  is  more  than  sufficient  to  determine  unambiguously  the 
spectroscopic  parameters. 

In  a  reaction  such  as  F  y  (a,  ny)Na  ,  the  sum  of 
spins  is  1  and  a  second  unquantized  parameter  is  thus  intro¬ 
duced  in  the  angular  correlation  formulae,  namely  the  ratio 
of  two  population  parameters  P(l)/P(0).  Since  three  para¬ 
meters  have  to  be  determined  now  from  experimental  data,  the 
analysis  may  not  lead  to  a  clear  cut  answer  in  some  cases. 


3 


Moreover,,  because  of  the  finite  size  of  the  neutron  counter, 
higher  magnetic  substates  can  be  populated  as  shown  later 
and  the  chances  of  getting  ambiguous  results  therefore 
increased . 


Thus,  before  performing  the  experiment,  it  had  to 
be  borne  in  minds  first,  that  no  theoretical  graph  of  the 
possible  angular  correlation  function  could  be  plotted 
because  of  these  two  unquantized  parameters;  consequently, 
it  would  not  be  easy  to  evaluate  the  validity  of  the  collected 
data;  secondly,  that  it  might  not  be  possible  to  limit  to 
one  spin  value  only  the  assignment  of  the  levels  to  be  studied. 

2 

Actually,  recent  computing  methods  based  on  "X 
fitting  programs",  (some  examples  of  which  are  given  by 
Pronko  (Pr  6 5)  )  allow  these  ambiguities  to  be  removed, 
and  it  was  known  that  as  a  last  resort,  it  would  be  possible 


to  have  recourse  to  them. 


. 
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1.2  TRIPLE  ANGULAR  CORRELATION 

In  the  past  few  years ,  the  theory  of  angular 
correlations  has  been  simplified  by  the  introduction  of  more 
powerful  mathematical  techniques  and  also  by  the  tabulation 
of  the  coefficients  involved  in  the  calculation  of  the 
theoretical  values. 

This  experiment  shows  one  example  of  triple  correla¬ 
tion  as  defined  by  Litherland  and  Ferguson  (Li  6l )  (Li  64) 
and  is  so  called  because  three  directions  are  involved  in 
the  correlation:  beam  axis  or  axis  of  quantization  (z-axis) 
and  directions  of  the  two  different  counters  with  respect 
to  the  beam  axis. 

After  recalling  briefly  the  general  formalism  of 
angular  correlation  theory.,  the  theoretical  formulae  of 
triple  angular  correlation  will  be  discussed  and  special 

IQ  22 

application  in  the  case  of  F  ~(a,  ny)  Na  reaction  will  be 
introduced. 


(a)  General  formalism  of  angular  correlations  for  aligned  nuclei. 

Considering  an  electromagnetic  radiation  produced 
by  deexcitation  of  an  axially  symmetric  state  with  angular 
momentum  a  and  emitted  in  the  direction  0  with  respect 
to  the  quantization  axis,  the  intensity  of  this  radiation  is 
given  by: 


* 
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W(9)  =  trace  ££  =  2  pk  e,  (2.1) 

a 

kn 

where  £  is  the  density  matrix  of  the  excited  level  and  its 
elements  are  related  to  the  population  of  the 
different  magnetic  sublevels. 

£  is  an  efficiency  matrix  so-called  because  it 
describes  the  actual  efficiency  of  detection  of 
the  radiation  in  the  direction  0. 

Because  the  excited  nucleus  is  considered  here  to 
be  aligned  (i.e.,  in  an  axially  symmetric  state)  one  can 
show  that  only  the  terms  n  =  0  are  non  vanishing  and  also 
that  the  density  matrix  is  diagonal. 

In  many  cases,  it  is  better  to  deal  with  tensors 
rather  than  with  matrices  because  tensors  transform  more 
simply  than  matrices  under  rotations.  These  tensors  are 
related  to  the  matrices  by  a  general  relationship  which, 
for  instance,  in  the  case  of  density  matrix  yields  the 
tensor  parameters i 

Pkn(aa,)=  2  (_i)aa  (a  a’  a-aT  |kn)<aa|p|a’a’> 


(2.2) 


6 

where  (a  af  a-aT  |  kn)  is  a  Clebsch-Gordon  coefficient. 

The  notation  used  here  is  the  notation  of  Devons  and  Gold- 
farb  (De  57) •  the  Greek  letters  represent  the  magnetic 
substates  of  the  state  tagged  by  the  equivalent  latin  letter. 

From  the  properties  of  the  Clebsch-Gordon  coeffi¬ 
cients  ,  it  can  be  seen  that: 

|  a  -  a*  |  k  ^  a  +■  a’ 

and 

n  =  a  -  a1 

Devons  and  Goldfarb  (De  57)  show  that,  because  of 
symmetry  conditions,  n  =  0  is  the  only  possible  value; 
thus  a  =  a’.  Moreover,  if  the  state  has  definite  parity  and 
if  no  polarization  is  present  in  the  incident  particle  and 
target  nucleus,  the  only  non  vanishing  tensor  parameters 
are  those  with  k  =  0  or  k  even.  The  eveness  of  k  implies 
that  the  density  matrix  elements  (population  parameters) 
are  equal  for  symmetrical  magnetic  sublevels;  that  is: 

P (a)  =  P( aT ) 


or 


<  aa  |  p  |  aTa*>  =  <a-a  |  p  |  a'-aT> 


(2.3) 
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For  each  particular  case,  the  angular  correlation 
formulae  may  be  derived  by  usings  for  instance,  the  proper¬ 
ties  of  tensor  products.  An  approximation  introduced  in 
the  evaluation  of  the  efficiency  tensor  parameters  is  to 
consider  a  plane  wave  propagating  along  the  z-axis  and 
detected  by  an  ideal  point  detector.  For  a  given  angle  0, 
the  rotation  of  the  efficiency  tensor  is  then  performed  and 
corrections  for  detectors  of  finite  dimensions  are  made. 

From  an  experimental  point  of  view  the  angular 
correlation  function  is  usually  defined  by : 

cl  Si 

W(e)  =  1  +  —  P„  (cos  e)  +  —  P,  (cos  0)  +  .  .  .  (2.4) 

ao  ao  4 

In  order  to  compute  the  different  a  coefficients 
by  the  method  proposed  by  Litherland  et  al  (Li  6l),  one  has 
to  calculate  different  coefficients  tabulated  by  Sharp  et  al 
(Sh  6l).  By  their  structure  these  coefficients  appear  as 
products  of  different  Clebsch-Gordon  and  Racah  coefficients 
and  a  concrete  example  will  be  given  in  the  case  of  the 
F1^  (any)  Na22  reaction. 

(b)  Triple  correlation. 

As  shown  in  the  diagram  hereafter,  a  particle  of 
orbital  angular  momentum  ^ ,  and  spin  s1,  impinging  on  a 


V 


8 

nucleus  of  spin  a  will  create,  when  a  reaction  occurs, 
an  intermediate  state  of  spin  b  which  will  then 
deexcite  by  emission  of  a  particle  of  orbital  angular 
momentum  and  spin  feeding  the  magnetic  sublevels  of 
a  state  c  which  decays  to  the  ground  state  by  emission  of 
an  electromagnetic  radiation  of  spins  L^,  (when  consider¬ 
ing  the  usual  case  of  two  multipoles  mixing) . 


s2 

-  c 

L1 ,  L2 

1 -  d 

The  counter  detecting  the  emitted  particle  is 
placed  at  0°  or  l80°  with  respect  to  the  z-axis  and  is 
first  assumed  to  be  equivalent  to  an  "ideal  point  detector. " 

Then,  strong  limitations  occur  in  the  number  of  magnetic 
sublevels  fed,  as  shown  by  Litherland  (Li  6l);  because  the  orbital 
angular  momenta  of  the  incoming  and  outgoing  particles  are 
perpendicular  to  the  z-axis,  only  the  intrinsic  spins  of 
the  various  particles  involved  are  contributing  in  population 
parameters . 


a 
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The  incoming  and  outgoing  channel  spins , 
respectively  and  S2  are  defined  by: 


S1  =  a  +  sx 


Sg  =  o  +  s„ 


(2.5) 


The  intermediate  state  although  it  has  no  sharp 
spin  is  determined  by: 


b  =  S-^  +  =  Sg  + 


(2.6) 


The  density  matrix  for  this  intermediate  state  b 
is  shown  to  be  (Li  6l)  : 


<bp  |  p  |  b'Pf>  =  -  2  (a  a  s^  |  S1  p)d  (Sn ,  p  0  |  b  p) 


ad 


T  K  1 


1 


x  (S1P^-[0  |  b  p)  <da  |  p  |  aa>  <s1  |  p  |  s1a1> 


11  r  “1  r 


<^10  |  p  |  l1  0> 


1 


x  <S1  &1  II  V  ||  b>  <S1l'1  ||  V  ||  b '  >  (2.7) 


If  the  states  are  assumed  to  be  unpolarized  then: 


.r 

10 


(2.8) 


(2.8-) 


And  by  considering  a  plane  wave  propagating  along  the 


z-axis  one  can  show  that: 


Ci-,0  |  p  |  i-jO  = 


1 


(2.9) 


The  "double  bar"  or  reduced  matrix  elements  are 
Introduced  by  the  perturbing  Hamiltonian  in  the  description 
of  the  absorption  causing  the  transition,  as  a  perturbation. 

From  the  properties  of  Clebsch  Gordon  coefficients, 
it  is  interesting  to  point  out  that: 


(3  =  a  +  ct 


(2.10) 


Similarly  the  density  matrix  of  the  final  state  c  is: 
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<C7S2^2^2 


From  this 
the  final 
involving 
spin  part 
the  third 


Now.,  from 


thus.,  the 


p  |  C7s2cr2i£0>  =  <07  |  p  |c7>  <s2a2  |  p  |  s2a2><i20  |  p  |i£0> 

2  (SppipO  |  bp)  (SpP^I 0  |  bp) 

bb'S2S2  ^  ^  ^ 


x  (c7S2a2  |  S2P)  (c7S2ct2  |  S2  p) 

x  <bp  |  p  |b'p>  x  <S2^2  ||  V  ||  b> 

<S2^2  ||  V  ||  b>* 


(2.11) 


equation.,  one  can  see  that  the  density  matrix  of 
state  is  the  product  of  three  density  matrices 
the  non  excited  final  state  and  the  orbital  and 
of  the  particle  exciting  this  state.  Also.,  from 
Clebsch  Gordon  coefficient.,  one  gets: 


P  =  7  +  °2  (2.12) 

(2.10)  and  (2.12) 

P  =  a  +  cr2  =  7  +  a2  (2.13) 

maximum  value  of  7  will  be  determined  by: 

7  =  max  (a  +  a  +  cj0)  =  a  +  s,  +  sQ  .  (2.14) 

'max  v  1  2 7  12 
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This  shows: 

1)  that  the  value  of  7  which  is  the  magnetic  quantum 
number  of  the  final  state  c  is  independent  of  the  spin  b 
of  the  intermediate  state  and  of  the  orbital  angular 

r  r 

momenta  i  2* 

2)  that  the  formula  for  the  angular  distribution  with 
respect  to  the  beam  of  the  electromagnetic  radiation  decaying 
from  a  state  of  spin  c  to  a  state  of  spin  d,  can  be  used 

for  angular  correlation,  i.e.,  when  the  intermediate  particle 

is  actually  detected  in  coincidence  with  the  electromagnetic 

radiation.  The  main  difference  between  the  two  ways  of 

investigation  comes  from  the  limitation  of  7  when  the 

max 

intermediate  particle  is  detected  at  0°  or  l80°. 

It  must  be  pointed  out  that  this  is  not  coinciden¬ 
tal,  since  basically  angular  distribution  formalism  and 
angular  correlation  formalism  are  similar.  Consequently, 
the  angular  distribution  or  angular  correlation  of  electro¬ 
magnetic  radiation  issued  from  a  state  c,  whose  population 
parameters,  P(y)  are 

P(y)  =  <cy  |  p  |  cy> 


and  leading  to  a  state  d,  is  given  by: 
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W(0)  =  2  (phase  vector) (density  matrix  of  state  c) 


x  (efficiency  matrix) 


(2.15) 


or 


Z-l  (L1cL2cJ>dk)QkPk(cos  0) 


(2.1 6) 


where 


f^-d  +  y+L-^  +  Lg  +  jj 


and  Lg  are  the  mixing  multipolarities  of 
the  electromagnetic  radiation 

5  is  the  multipole  mixing  ratio  (see  remark 
number  (ii)  hereafter) 

Z1  coefficient  tabulated  by  Sharp  et  al  (Sh  6l) 

IQ  22 

p  =  0,  1  or  2  (see  case  of  F  ^(any)Na  ) 

is  the  attenuation  factor  due  to  the  fact 
that  the  radiation  is  actually  detected 
by  finite  size  detectors.  These  values 
have  been  tabulated  for  different  sizes 
of  detectors  and  depend  on  the  energy 
of  the  7  ray  and  the  crystal -source  distance. 
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Remarks : 

i)  In  the  case  of  angular  correlation  with  detection 
of  outgoing  particles  at  0°  or  l80°,  it  has  been  assumed 
that  the  detector  is  ideal  (i.e.,  very  small);  actually, 
there  is  a  finite  size  of  detector  which  cannot  be  neglected. 
Litherland  and  Ferguson  (Li  6l )  have  shown  that  if  the  half 
angle  subtended  by  the  counter  is  £  (in  radians),  the 
substates  beyond  7max  are  actually  populated  and  the  sub¬ 
states  (y  +  n)  are  populated  with  a  population  proportional 
max 

2n  ■ 

to  i  .  These  calculations  suppose  that  the  detector  is 
on  the  axis  of  the  beam.  When  this  is  not  the  case,  i.e. 
detector  at  0°  but  its  axis  being  perpendicular  to  the  beam 
direction,  it  is  reasonable  to  assume  that  these  results 
are  still  valid  although  the  definition  of  the  angle  I*  is 
less  easy. 

ii)  If  the  angular  correlation  is  being  expressed  as: 

w(0)  =  WLlL4e)  +  26  WLlL2  +  62  WL2L2(e)  L1  <  L2 

then  6  is  defined  by: 

where  <d  ||  II  c>  and 
<d  |  |L2  ||  c>  are  the  reduced 
matrix  elements  of  the 
multipole  radiation. 


amplitude  of  <d  ||L2  ||  c> 

amplitude  of  ^  =  <d  ||Li  ||c> 


. 
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111)  Since  angular  distribution  and  angular  correlation 

are  identical,  an  alternate  method  can  be  used  in  the  case 

where  neutrons  are  the  outgoing  particles.  Effectively, 

when  neutrons  have  a  low  energy,  or  in  other  words,  when 

the  state  fed  by  neutrons  is  just  above  threshold,  a 

limitation  of  y  is  again  possible  since  low  energy 

nicix 

neutrons  have  zero  or  low  angular  momentum.  However,  this 
method  is  only  possible  when  a  resonance  occurs  just  above 
threshold . 


IQ  22 

(c)  Application  to  F  ^  (a,ny)Na  reaction. 

Due  to  the  limitations  brought  by  recent  experi¬ 
ments  (see  section  1.3  c),  the  calculations  for  the  unknown 
spins  of  the  2  states  892  kev  and  1.532  Mev  have  been 
carried  out  assuming  all  possible  values  for  the  unknown 
spins  between  1  and  4. 

A  very  unlikely  0  spin  because  of  high  order  mul¬ 
tipolarity  radiation  would  yield  both  angular  distribution 
and  angular  correlation  isotropic. 

In  the  case  of  spin  1  the  diagram  of  the  reaction 


is  the  following: 


Since  the  parity  is  known.,  one  should  expect  in 
this  case  a  predominant  E2  radiation  with  weak  admixture  of 
M3.  As  shown  in  section  1.3  c,  the  Weisskopf  estimates  give 

I  g  fs 

a  ratio  E^/M^  varying  from  10  to  10  for  the  2  states 
studied. 


Using  the  formula  (2.l6),  the  calculations  are 
then  carried  out  for  P(0)  and  P(l)  since  0  and  ±1  are  the 
only  magnetic  sublevels  which  can  be  populated. 

The  complete  expression  of  W(0)  results  as  follows: 


w(0)  1  +  0.144  (P(O)-P(l)}  +  1.170  (U0)-P(1))5-0.250(P(0)-P(1)}5 


(l+o  )  (P(0)  +  2P ( 1 ) ) 


X  P„  (cos  0) 


w(e) 


0.144  +  1.1705  -  0.25052  [p(0)  _  p 

(1  +  62) (P(0)  +  2P ( 1 ) ) 


(cos  d) 


or 


. 
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If  the  population  parameters  are  equals  one  can 
check  that  equally  populated  substates  (or  in  other  words, 
non  aligned  nuclei)  will  be  revealed  by  an  isotropic 
distribution  W (0)  =  constant. 

In  order  to  make  the  interpretation  of  experimental 

results  easier,  the  graph  &2/&0  versus  ^  for  different  values 
2  P  (l ) 

of  x  defined  as  p is  plotted. 

This  has  been  done  for  all  the  four  possibilities 
and  assuming  s  -  wave  neutrons  (see  appendix  Al).  As  may 
be  expected,  some  cases  would  lead  to  ambiguous  results 
but  in  the  most  pessimistic  case,  the  number  of  possible 
assignments  should  be  reduced  from  5  to  2,  for  instance  a 
negative  value  of  a.^/a.Q  would  limit  the  possible  spins  to 
2  and  3>  a  positive  value  to  2  and  4,  etc.  .  .  . 
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1.3  The  Na  nucleus. 

a.  Theoretical  aspects. 

Many  models  are  available  today,  which  try  to 
describe  nuclear  structure.  They  are  more  or  less  successful 
in  explaining  the  physical  properties  of  nuclei  and, 
although  even-even  nuclei  and  odd-even  nuclei  can  often  be 
understood  theoretically,  it  appears  that  no  theory  ade¬ 
quately  explains  the  physical  properties  of  all  odd-odd 
nuclei.  In  other  words,  closed  shell  nuclei  or  nuclei  which 
can  be  described  by  a  one  particle  state  are  the  easiest 
to  explain  theoretically  whereas  nuclei  such  as  odd-odd 
nuclei  whose  configuration  is  far  from  closed  shells  cannot 
be  so  easily  comprehended. 

Two  main  approaches  have  been  taken  into  considera¬ 
tion  by  the  theoreticians.  The  first  one  was  logically 
introduced  by  the  previous  study  of  atomic  structure : 
since  atomic  phenomena  were  satisfactorily  explained  by 
considering  each  electron  separately,  one  tried  to  explain 
the  nucleus  considering  it  as  an  aggregate  of  nucleons 
obeying  similar  laws  of  quantization.  This  gave  rise  to  the 
shell  model  which  leads  to  a  satisfactory  matching  with  many 
experimental  results,  although  it  fails  to  account  for  the 
large  nuclear  quadrupole  moments  and  spheroidal  shape  which 
many  nuclei  have.  The  second  approach  then,  was  to  consider 


■ 
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the  nucleus  in  two  parts;  a  core  whose  shape  would  keep 
almost  constant  in  time  and  some  extracore  nucleons.  This 
second  model  is  usually  referred  as  the  "collective  model" 
although  the  terminology  is  not  clear. 

Finally,  a  kind  of  compromise  was  attained  with 
the  work  of  Nilsson  and  the  Danish  school:  the  Nilsson 
model,  referred  sometimes  as  the  "unified  model"  could  be 
defined  as  a  more  sophisticated  shell  model  in  which  collec¬ 
tive  effects  are  taken  into  account.  The  energies  of  the 
single  particle  levels  are  function  of  the  deformation 
parameter  of  the  core.  A  comprehensive  series  of  papers 
can  be  found  in  the  Danish  publications  from  1950  to  i960, 
the  most  important  of  them  being  probably  (Ni  55) • 

22 

From  the  viewpoint  of  the  shell  model,  Na  has 
the  configuration  p(l  d^2)J  n(l  ^5/2’  3  °"^her  wor(is, 

outside  the  closed  p^-proton  shell,  there  are  3  protons  in 
the  dj-/2  shell  and  in  order  to  complete  this  1  d^g  shell, 
one  would  need  to  add  3  more  protons.  The  situation  is 
identical  for  the  neutron  shells.  This  also  means  that  if 
the  last  filled  shell  limits  an  inert  core,  one  would  have 
to  solve  a  six  body  problem  to  get  the  state  function  of  the 
nucleus.'  However,  under  some  simplifying  assumptions,  the 
shell  model  allows  a  satisfactory  determination  of  some 
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physical  constants  like  nuclear  magnetic  moments  and  ground 

state  spin  assignments  even  in  such  an  extreme  case.  The 

Schmidt  prediction  for  nuclear  magnetic  moments,  although  it 

deals  with  an  extreme  single  particle  model,  can  be  applied 

22 

to  odd-odd  nuclei  and  in  the  case  of  Na  yields  a  value  of 

1.73  nuclear  magnetons  (Ra  57)  which  compares  closely  to 

the  known  experimental  result:  1-75  nuclear  magnetons.  In 

22 

order  to  determine  the  ground  state  spin  of  Na  ,  Nordheim's 
rule  as  modified  by  Brennan  and  Bernstein  (Br  60)  can  be 
applied.  As  shown  by  Brennan  and  Bernstein,  the "weak  rule" 
of  Nordheim  can  be  replaced  by  a  stronger  rule  (R2).  If 
J1  and  (or  J 2  and  are  the  single  particle  total  and 

orbital  angular  momenta  obtained  from  the  adjacent  odd  A  nuclei 
using  standard  shell  model  assignment,  then  the  ground  state 
spin  J  is  given  by: 


J  =  I  J1  ±  J2  I  for  Ji  =  ^1  ±  i 

and  J  ^  ±  i 

pi  PS  3 

Since  Ne  and  Na  both  have  ground  state  spin  ^  +  ,  this 

leads  to  two  possible  values  for  the  ground  state  spin  of 
22 

Na  ,  namely  3  and  0. 
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Before  introducing  the  collective  aspects  of  the 
nucleus,  it  may  prove  useful  to  point  out  that  the  formalism 
used  hereafter  is  that  introduced  by  the  Danish  school  and 
can  be  recalled  briefly  by  the  following  diagram: 


Fig.  1.3.1 


where  Z 

7 

7 

-f 

J 


is  the  quantization  axis  and  ZT  an  arbitrary  axis 
is  the  total  angular  momentum  of  the  nucleus,  M  and 

K  its  respective  projections  on  the  Z  and  ZT  axes, 
is  the  total  angular  momentum  of  the  core 
is  the  total  angular  momentum  of  the  external 

particle (s)  and  Q  the  projection  of  J  on  the 


Z ’  axi s 
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In  the  case  of  axially  symmetric  nuclei ,  there  is 
no  contribution  of  the  core  angular  momentum  along  ZT  if 
Z T  is  symmetry  axis;  therefore  K  =  £2.  This  can  be  shown  from 
classical  mechanics  considerations  (see  e.g.  (Pr  63)  p.  237 
and  256) . 


As  seen  earlier,  both  unpaired  proton  and  neutron 

22 

of  Na  are  in  the  the  shell  model  representa¬ 

tion.  Nilsson  (op.  cit.)  has  shown  that,  by  considering  a 
non  spherical  shell  potential,  this  d^2  level  splits  into 
3  "orbits"  as  the  deformation  increases.  These  orbits 
corresponding  to  single  particle  states  are  labelled  by  the 
Nilsson  number,  the  intrinsic  spin  value  £2  and  the  parity 

of  the  state.  According  to  this  model,  one  should  expect 

22 

for  the  ground  state  of  Na  both  odd  neutron  and  proton  in 
Nilsson  orbit  7(j|+)  j  (Ar  6 5)  j  hence  two  possible  spin 
assignments  corresponding  to  parallel  and  antiparallel 
coupling  of  the  intrinsic  spins: 


and 


£2  +  £2 


p  n 


J77  =  l  £2  -  £2  |  =  0+  . 

f  p  n  1 


It  can  be  shown  that  parallel  coupling  is  the  most  probable 
(Na  65)  and  this  fact  is  confirmed  by  experiment  as  shown  in 

table  1  of  a  paper  by  Rakavy  (Ra  57).  Experimentally,  the 

22  + 
ground  state  spin  of  Na  is  known  to  be  3  . 


' 
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Rakavy  (op.  cit.)  has  shown  that  Nilsson  model 
could  be  more  extensively  applied  in  the  region  between 
A  =  18  and  A  =  28.  According  to  his  calculations  the  defor¬ 
mation  in  the  nuclei  following  0^  increases  very  rapidly 
and  for  A  =  18  already  attains  a  value  characteristic  of 
heavy  nuclei  revealing  rotational  spectra*  i.e.,  for  nuclei 
whose  A  varies  from  l6o  to  180.  The  calculated  deformation 

(which  is  known  as  6  or  e  in  the  literature)  reaches  a 

20 

maximum  for  Ne  then  stays  almost  constant  until  A  =  25  is 
reached*  then  decreases.  Evidence  of  the  fact  that  rotational 
spectra  are  found  in  the  region  A  =  18  to  A  =  28  is  shown  in 
a  summary  of  pertinent  work  (Gov  6o) .  Mottelson  and  Nilsson 
(Mo  59)  pointed  out  that  the  deformation  parameter  6 

decreases  with  increasing  A  because*  in  the  model  which  is 

_1 

studied,  the  deformations  are  fundamentally  of  order  A  3  . 

Commenting  upon  Rakavy' s  paper*  Mottelson  and  Nilsson  mention 

that  the  experimental  data  are  not  sufficiently  complete  to 

provide  a  very  critical  test  of  the  coupling  scheme  as 

interpreted  by  Rakavy*  although  it  does  provide  a  satisfactory 

basis  for  interpreting  most  of  the  available  evidence  in 

this  region  (A  =  18  to  A  =  28).  Finally*  it  is  interesting 

to  point  out  that  Varshalovich  (Va  64)  gives  some  evidence 

22 

to  support  the  fact  that  Na  is  probably  strongly  deformed* 
by  studying  the  first  two  excited  states  of  this  nucleus. 
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As  will  be  shown  hereafter,  only  3  states  of  Na 
have  been  assigned  a  definite  spin  value  and  it  is  therefore 
difficult  to  compare  this  nucleus  to  any  model.  However, 
for  the  reasons  explained  earlier,  the  rotational  aspect 
of  the  nucleus  seems  to  be  the  most  applicable  and  will 
thus  be  taken  into  consideration. 


22 

b.  A' review  of  previous  work  on  Na 

At  the  present  time,  23  excited  states  have  been 
22 

found  for  Na  (Ar  65)*  but  only  three  of  them  have  been 
assigned  a  definite  spin  and  parity  whereas  definite  parity 
and  probable  values  for  spin  have  been  assigned  for  the  next 
two.  (fig.  1.3.2. ) 

MeV 


(l+)4  4+ _ 

T=1,0+  j+  _ A 
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Fig.  1.3.2.  Diagram  illustrating  the  low-lving  energy  levels 

and  gamma-ray  branching  in  Na22  after  Arnell  and 
Wernbom-Selin  (Ar  6 5) • 
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By  studying  the  pick-up  reaction  Na  J  (d,t)  Na  , 
Vogelsang  and  McGruer  (Vo  58)  showed  that  the  first  three 
T  =  0  excited  levels  have  positive  parity  and  that  their 
spin  value  is  limited  to  4. 

Temmer  and  Heydenberg  (Te  58)  found  the  first 
T  =  1  state  at  660  kev.  This  state  is  a  0+  and  decays  to 
the  ground  via  the  587  kev  state. 

oh  22 

Pellegrini  (Pe  6l )  studied  the  Mg  (d,a)  Na 
reaction  and  tried  to  fit  the  angular  distributions  of 
alpha  particles  with  D.W.B.A.  theory  and  Butler  theory. 

For  the  1.532  Mev  State  the  best  fits  were  obtained  with 
Butler  calculations  for  i  -  0  and  $  =  4,  thus  limiting  the 
possible  values  of  J 71  to  1+  or  3+ ,  4+,  5+*  For  the  892  kev 
state,  the  best  alpha  particle  distributions  fit  was 
obtained  with  D.W.B.A.  calculations  for  A  =  0.  Since  this 
fit  was  poor,  Pellegrini  did  not  give  a  possible  assignment 
for  this  state.  However  Endt  and  Van  der  Leun  (En  62) 
referred  to  Pellegrini's  work  and  concluded  that  a  probable 
spin  assignment  for  the  892  kev  state  was  (l+). 

22 

In  the  past  few  years,  research  on  Na  was  oriented 
much  more  towards  the  discovery  of  the  second  T  =  1  state.  As 
seen  previously,  the  first  T  =  1  state  is  the  660  kev  level 
which  is  the  analogue  of  the  ground  state  of  Ne2^  (T  =  1,  J7r=0+). 
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Since  the  first  excited  state  of  Ne  is  the  1.28  Mev  level, 

22 

the  second  T  =  1  state  of  Na  is  expected  at  an  excitation 
energy  of  about  1.94  Mev.  None  of  the  three  levels  at  1.942 
1.95  and  1.988  MeV  exhibits  characteristic  properties  of  a 
T  =  1  state  and  the  problem  remains  unsolved. 


c.  Predictions. 
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As  seen  earlier,  Na  can  be  considered  as  a 
strongly  deformed  nucleus  exhibiting  a  pure  rotational 
spectrum.  Rakavy  (Ra  57)  gives  0.47  as  the  calculated 
equilibrium  deformation  and  predicts  a  value  of  1.78  n.m. 
for  the  nuclear  magnetic  moment,  known  to  be  1.75  n.m. 
experimentally.  Considering  that  there  are  2  odd  particles 

o 

in  the  same  Nilsson  orbit  7(75+)  j  three  low  lying  bands  can 
be  expected  according  to  Rakavy: 


K  =  2Q  =  3+ 

with 

T  = 

0 

(the  lowest  level 
this  band  is  the 
ground  state) 

w 

II 

0 

+ 

with 

T  = 

0 

(  Jlr=l+,3+,5+...) 

II 

0 

+ 

with 

T  = 

1 

(  J7r=0+,2+,4+...) 

Paul  (Pa  65)  proposes  a  more  complete  level  scheme 
by  considering  all  the  possible  transitions  for  the  two  odd 
particles  to  higher  Nilsson  levels.  He  thus  gets  the  differ 
ent  following  combinations: 
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1)  both  proton  and  neutron  in  Nilsson  orbit  7(|+).  This 
gives  bands  of  K  =  3+  and  0+ . 

2)  one  particle  in  level  9(75+)  and  the  other  in  level 
7(^+),  this  gives  four  bands:  two  of  K=l+  and  two  of  K=2+. 

3)  one  particle  in  level  5(^t)  and  the  other  in  level  7(^+) 
two  bands  of  K=4+  and  two  of  K=l+. 

Once  more*  it  is  necessary  to  point  out  that  these 
predictions  are  rather  speculative  since  the  experimental 
investigations  have  not  given  enough  information  so  far. 

It  has  been  seen  that  all  the  possible  spin 
assignments  for  the  states  892  kev  and  1.532  MeV  were  limi¬ 
ted  to  a  maximum  value  of  4  with  positive  parity  and  that  a 
probable  spin  assignment  for  the  892  kev  state  would  be  1+. 

If  one  considers  the  simple  case  of  single  particle  transition, 
these  possibilities  lead  to  two  different  multipole  radiation 
mixtures  Ml  and  E2  for  a  transition  2+  (3+  or  4+)  to  3+j 
Ml  being  predominant,  E2  and  M3  for  a  transition  1+  to  3+ 

E2  being  predominant.  As  pointed  out  by  many  authors  (see 
(Pr  63)  p.  341  for  instance),  the  Weisskopf  estimates  for 
single  particle  transitions  may  be  off  by  a  few  orders  of 
magnitude.  However,  these  estimates  generally  give  a  basis 
for  discussion  and  it  may  prove  useful  to  calculate  them. 
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For  single  particle  transition,  the  mixing  ratios 
of  multipole  radiation  are  given  by: 

Ml  2.8  in5  1  1 

E2  1U  ^V3  e2 

7 

for  transition  2+,  3+j  to  3+ ,  and 

E2  1,8  n  J  1  1 

M3  5  7+73  f2 

7 

for  transition  1+  to  3+ > 

where  A  is  the  atomic  number  and  E^  the  energy  of  the 
transition  gamma-ray  in  MeV. 


For  the  1.532  MeV  transition,  one  gets: 


Ml 

E2 

E2 

M3 


1260 

6.5  x  105 


For  the  .892  MeV  transition,  one  gets: 


Ml 

E2 

E2 

M3 


3500 

1.8  x  10^ 
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It  is  interesting  to  point  out  that,  according  to 
Morpurgo  (Mo  63)  *  El  transitions  are  forbidden  for  the  AT  =  0 
transitions  in  self  conjugate  nuclei:  N  =  Z  or  =  0. 

Under  the  same  conditions.  Ml  transitions  are  expected  to 
be  reduced  by  a  factor  100  to  150.  This  rule  does  not 
affect  E2  transitions.  Experiments  have  shown,  however, 
that  electric  quadrupole  transitions  are  enhanced  compared 
to  the  single  particle  approximation  when  the  nucleus 
exhibits  collective  effects.  (see  (Mo  63)  e.g.). 

Thus  in  the  case  of  transitions  2+  (3+  or  4+)  to 

3+,  one  can  expect  strong  Ml  E2  mixing  in  self  conjugate 

nuclei.  As  shown  by  Morpurgo  (op.  cit.),  the  quadrupole 

transition  may  even  become  predominant  in  the  case  where 

E  >  16/A2/3  (in  MeV);  this  would  correspond  to  E^  >  2  MeV 
22 

for  Na  . 
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CHAPTER  II.  EXPERIMENTAL  ASPECTS 

2.1  Experimental  apparatus, 
a.  Correlation  table. 

The  beam  of  alpha  particles  was  accelerated  by  the 

■jf 

5.5  MV  electrostatic  accelerator  of  the  Nuclear  Research 
Centre.  The  correlation  table  was  the  same  as  used  by 
J.  G.  Pronko  (Pr  65)  but  some  changes  have  been  introduced 
in  the  beam  path  (fig.  2.1)  and  the  collimator  has  been 
adjusted  to  fit  exactly  the  inside  diameter  of  the  tube. 

The  targets  were  prepared  at  the  Nuclear  Research 
Centre  laboratory  and  consisted  of  a  0.002  inch  thick 
platinum  backing  on  which  CaF^  had  been  evaporated.  In  the 
first  series  of  runs,  a  gold  backing  was  added  for  the  pur¬ 
pose  of  better  heat  transfer.  For  reasons  explained  later, 
this  gold  backing  was  removed  in  the  second  series  of  runs. 

The  target  was  cooled  by  compressed  air  stream  during  the 
experiments . 

Model  CN  Van  de  Graaff,  High  Voltage  Engineering  Corporation, 
Burlington,  Massachusetts,  U.  S.  A. 
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CM 


cn 


An  illustration  of  the  beam  path. 
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Neutrons  were  detected  by  a  pulse  shape  discrimi- 
nator  probe  .  The  1-J-"  long  x  1^"  diameter  liquid  scintillator 
is  mounted  on  a  Philips  56  AVP  photomultiplier. 

Different  sodium  iodide  crystal  sizes  (2M  x  2", 

3"  x  3"  and  4"  x  4")  and  different  scintillator  photo- 


multiplier  combinations 

were 

used 

for  the  y-detection: 

4”  x  4" 

Nal 

(Tl) 

with  RCA  7046 

4"  x  4" 

Nal 

(Tl) 

with  58  AVP  (XP  1040)  Philips 

3"  x  3” 

Nal 

(Tl) 

with  54  AVP  Philips ** 

A  schematic  diagram  of  the  electronic  equipment  is 
given  in  illustration  (2.2).  This  is  equivalent  to  a  fast- 
slow  coincidence  unit  which  can  be  analyzed  as  follows: 

Fast  side  of  the  coincidence  unit: 

From  the  last  dynode  of  the  y-detector,  the 
positive  pulse  (2  volts,  60  microseconds)  is  limited  then 

*)f  -X"  ■Jf 

amplified  by  a  wide  band  distributed  amplifier  before  it 
enters  the  time  converter. 

"P.S.D."  Model  NE  5553*  Nuclear  Enterprises  Ltd.,  Winnipeg. 

Model  "scintibloc" .  Quartz  et  silice,  8  rue  dT Anjou, 

Paris,  FRANCE. 

■)f 

Model  460  AR,  Hewlett  Packard  Corporation,  Palo  Alto, 
California,  U.  S.  A. 


,•*0'  . 


33 


Inve 

;rt . 

H 

— i 

to 

C t 

< 

E 

TO 

Pulse 

Snaper 

*-  • 

to  (_> 

A. 

TO 

u 

c 

"D 

TO 

3 

4-1 

»- 

•r-< 

TO 

TO 

r— < 

> 

E 

C- 

C 

•*“4 

E 

o 

H 

< 

CJ 

Q 

.to 

r 

tu 


u 

X.  H 

OJ~ 

O  TO 

o  ^ 

co 

o 

ao 

^ — 

s 

O 

OS 

OS 

OJ 

H 

H 

u 

PS 

tJ 

o 

a- 

to 

to 

H 

CO 

TO 

r— 4 

JO 

>> 

TO 

TO 

•r*4 

r— 4 

u 

TO 

TO 

Q 

> 

TO 

»— 1 

JO 

TO 

TO 

•H 

—H 

o 

TO 

TO 

Q 

> 

TJOOI 


CM 

CM 

cn 


to 


CM 

to 


An  illustration  of  the  electronics. 
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From  the  dynode  14  (fast  output)  of  the  neutron 
detector,,  the  positive  pulse  is  amplified  by  a  double  delay 
line  amplifier  and  limited  by  a  clipping  circuit. 

The  two  pulses  (n  and  y)  entering  the  time  to 
amplitude  converter  (Ne  60)  are  +5  volts  in  amplitude  and 
150  nanoseconds  long.  The  output  of  the  time  converter  is 
a  negative  pulse  of  4  volts  maximum  amplitude  with  a 
risetime  of  approximately  100  nanoseconds. 

This  fast  coincidence  which  is  used  as  time  base., 

is  then  fed  into  a  single  channel  analyzer  whose  purpose  is 

to  determine  the  resolving  time  of  the  coincidence  unit  by 

setting  a  window  on  the  time  converted  signals.  This  window 

was  set  in  order  to  accept  only  pulses  whose  amplitude  was 

greater  than  3  volts.  Knowing  that  4  volts  corresponds  to 

a  time  of  150  nanoseconds*  this  fixes  the  resolving  time  at 

approximately  40  nanoseconds.  The  window  was  regularly 

checked  and  no  apparent  change  in  its  width  was  noticed. 

A  shaper  following  the  single  channel  analyzer  generates  an 

output  pulse  1  microsecond  wide  and  1  volt  high.  This  pulse 

is  put  into  coincidence  with  the  neutron  pulse  from  the 

pulse  shape  discriminator  output  and  the  coincident  signal 

provides  the  gate  pulses  for  the  1024  channel*  two  parameter* 

* 

pulse  height  analyzer. 

Model  ON  1024*  Technical  Measurement  Corporation*  44l 
Washington  Avenue*  North  Haven*  Connecticut*  U.S.A. 
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Slow  side  of  the  coincidence  unit: 

The  "slow"  pulses  from  dynode  10  for  the 

gamma-detector  and  for  the  linear  output  for  the  neutron 
detector  are  fed  into  separate  lines  including  linear 
preamplifier  (GO  60)  and  double  delay  line  amplifier  .  After 
amplification,  the  pulses  traverse  a  variable  delay  box  be¬ 
fore  entering  their  respective  inputs  at  the  pulse  height 
analyzer. 

Delays  are  set  in  the  different  lines  to  cancel 

the  difference  in  cable  lengths  and  to  match  the  arrival  of 

slow  and  fast  signals  at  the  pulse  height  analyzer.  In 

22 

order  to  set  these  delays,  a  Na  source  is  placed  between 
the  two  detectors  at  180°  one  from  the  other  and  the  simul¬ 
taneous  emission  of  the  two  annihilation  gamma-rays  is  taken 
as  time  reference. 

Singles  spectra  can  be  recorded  by  switching  the 
input  selectors  of  the  pulse  height  analyzer  to  anticoinci¬ 
dence,  thus  ignoring  the  fast  gating  pulses. 

Transistor  double  delay  line  amplifier  AEP  1436  Atomic 
Energy  of  Canada  Limited,  Chalk  River,  Ontario  (cf.  GO  6o) . 

Type  2011,  AD-YTJ  Electronics  Limited,  Passaic,  New  Jersey, 

U.  S.  A. 


36 


For  monitoring  purposes,,  the  PSD  output  correspond¬ 
ing  to  neutron  pulses  is  connected  to  a  scaler.  A  gamma-ray 
monitor  was  also  used  during  some  runs.  It  consists  of  a 
2"  x  2"  crystal  and  its  associated  photomultiplier,  a  preampli¬ 
fier,  a  double  delay  line  amplifier  and  a  double  channel 
analyzer  before  the  scaler.  The  purpose  of  this  second  monitor 
was  first  to  check  the  monitoring  by  neutrons  and  also  to 
have  a  reliable  normalization  when  the  neutron  detector  was 
moved  in  order  to  determine  its  effect  on  the  gamma-detector 
for  small  angles.  (See  "corrections"  (section  2.4)  ). 


2.2  Preliminary  experimental  results. 

A  test  run  was  performed  a  few  months  before  the 

1  Q 

experiment  started.  A  previous  study  of  F  +  a  reaction 
by  Freeman  and  Mani  (Fr  64)  has  shown  that  the  1.532  MeV 
state  of  Na  is  excited  at  a  resonance  of  4.7  MeV  bombarding 
energy.  One  of  the  spectra  obtained  during  this  test  run  is 
shown  in  fig.  2.3  and  corresponds  to  a  bombarding  energy  of 
approximately  4.7  MeV.  As  can  be  seen,  three  competing 

* 

reactions  are  occurring:  inelastic  scattering  F1^  (a,aT)F^  , 

22 *  * 

(a,p)  reaction  leading  to  Ne  (Q  =  1.675  MeV)  and  finally 

(a,n)  reaction  leading  to  Na^  (Q  =  -  1.9^2  MeV)  . 


587  keV 

Ca.n' 
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50  100  Channel  number. 

Fig.  2.3:  Ungated  spectrum  for  a  bombarding  energy  of  4.7  MeV . 
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When  the  experiment  actually  started,  it  was  noticed 
that  both  gamma-ray  and  neutron  yields  were  high  but  apparent¬ 
ly  there  was  no  overloading  of  the  electronics  for  a  beam 

current^  0.5  m-A.  However,  the  target  which  had  been  used 

o 

for  the  test  run  (  nj  100  pg/cm  of  CaFg  on  a  0*002  inch  thick 
platinum  backing)  appeared  to  burn  out  easily  with  such  a 
current,  although  it  was  cooled  by  air.  It  was  then  decided 
to  use  a  thicker  target  (aJ  300  (ig/cm  of  CaF^,  same  backing) 
and  to  mount  this  target  with  a  1/32  inch  thick  gold  backing 
for  better  heat  transfer.  The  beam  current  on  the  target 
was  kept  at  a  value  of  about  0.5  and  the  corresponding 
yield  gave  approximately  5*500  counts  per  second  in  the 
gamma-detector  and  approximately  5*000  counts  per  second  in 
the  neutron  detector,  both  values  supposed  to  be  acceptable. 
Since  the  target  then  was  about  15  times  thicker  than  the 
target  used  by  Freeman  and  Mani,  the  sharp  resonances  they 
had  observed  could  not  be  expected.  Actually,  the  resonances 
observed  with  the  300  |rg/cm  target  were  about  200  kev  wide. 
After  a  series  of  runs  in  the  region  4.5  MeV  to  5.5  MeV, 
the  energy  of  4.79  MeV  (from  nuclear  magnetic  resonance 
field  reading)  was  chosen  because  it  produced  maximum  excita¬ 
tion  of  the  1.532  MeV  state, 


. 
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The  main  problem  which  had  to  be  faced  then  was  a 
large  gain  shift  attaining  maximum  values  of  more  than  30$ 
per  day.  By  feeding  a  pulse  of  constant  amplitude  into  the 
gamma  preamplifier,  it  was  found  that  it  was  not  the  elec¬ 
tronics  but  the  detector  itself  which  was  faulty.  Different 
detectors  and  different  crystal-photomultiplier  combinations 
were  tried  but  they  showed  the  same  defect.  However,  since 
the  gain  shift  was  noticeably  smaller  with  a  4"  diameter  x 
4  "  long  crystal  mounted  on  an  XP  1040  photomultiplier,  an 
angular  correlation  experiment  was  tried.  It  appeared  that, 
although  the  coincidence  spectra  were  well  defined,  the 
results  were  not  reliable  after  all  corrections  had  been 
applied.  Moreover,  these  runs  would  only  concern  the  587  kev 
and  892  kev  states  since  the  statistics  were  too  poor  to 
investigate  the  1 , 532  MeV  state. 

As  shown  in  fig,  2.4,  it  had  been  noticed  that  for 
small  angles  (0°  to  30°)  two  unidentified  peaks  appeared 
which  corresponded  to  approximately  220  and  3^0  kev,  The 
second  peak  at  340  kev  disappeared  almost  completely  (fig,  2,5) 
when  a  2  inch  thick  piece  of  lithium  doped  paraffin  was 
placed  in  front  of  the  gamma  detector,  showing  that  this 
secondary  effect  was  due  to  neutrons,  Also,  the  sodium 
iodide  crystal  waB  getting  activated  by  the  high  neutron  yield. 


220  keV 


4o 


uoDzi-y) 


CLUiot; 


U  I  <  Z  Z  UJ  — I 


50  100  150  200  250  30C 

Channel  number. 

Fig.  2.4:  Neutron  gated  spectrum  for  a  bombarding  energy  of  5.01  MeV. 

(4"  x  4"  crystal  at  30  degrees  from  the  beam  axis) 
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After  a  few  hours  of  running,  the  resolution  of  the  detector 
was  much  poorer  than  initially.  Then,  with  no  beam  on  the 
target,  a  background  spectrum  was  taken  which  showed  high 
activation.  Resuming  the  experiment  after  a  few  minutes 
with  no  beam  on  the  target,  it  was  found  that  the  resolution 
of  the  detector  had  improved  drastically. 

Finally,  a  series  of  tests  with  radioactive  sources 
showed  that  the  counting  rate  of  the  3"  x  3"  "scintibloc" 
had  to  be  limited  to  a  maximum  value  of  3*000  counts  per 
second  in  order  to  limit  the  gain  shift  to  less  than  2 % 
per  day. 


2.3  Experimental  procedure. 

a)  In  order  to  improve  the  experimental  results,  the 
following  conditions  were  established: 

i)  limitation  of  the  counting  rate  to  3  x  10  counts 
per  second  for  the  gamma  detector. 

ii)  low  beam  current  limited  to  a  maximum  value  of 

0.2  (jlA . 

iii )  use  of  an  average  thin  target,  approximately  100 
|i.g/cm2  of  CaF2  on  a  0.002"  thick  platinum  backing.  No  gold 
heat  conducting  back  was  used. 
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A  2.5"  thick  lithium  doped  paraffin  cap  was  placed  in  front 
of  the  gamma  detector  to  minimize  the  activation  due  to 
neutrons  and  also  fulfill  the  condition  (i)  limiting  the 
gamma  counting  rate.  It  was  also  decided  to  take  one 
singles  spectrum  at  90°  after  each  run  for  checking  and 
monitoring  purposes. 


b)  Kinematics  of  the  reaction. 

The  kinematic  variables  of  the  reaction  were 
computed  using  a  program  written  for  the  LGP  30  computer  of 
the  Nuclear  Research  Centre.  Prom  the  results  of  the 
computations j  the  energy  threshold  for  each  excited  state 
and  the  energy  of  neutrons  emitted  at  0°  for  a  given  bom¬ 
barding  energy  were  obtained.  The  Q-value  of  the  reaction  is 
-1.949  MeV  and  the  nuclear  masses  are  as  follows: 

M-,  =  4.003876  amu  1 

>  input  channel 

M2  =19.00444  amu  J 


=  I.OO8986  amu 

M^  =  22.00143  amu 


output  channel 


1  amu  =  931 • l4l  MeV 


This  resulted  in  the  following  table: 


44 


Level  in  Na^ 

Threshold 

Energy 

2° 

Level  in  Na 

Threshold 

Energy 

ground  state 

587  kev 

660  kev 

892  kev 

2.359  MeV 

3.070  MeV 

3.159  MeV 
3.439  MeV 

1.532  MeV 
1.942  MeV 

1.950  MeV 
1.988  MeV 

4.214  MeV 

4.711  MeV 
4.720  MeV 
4.766  MeV 

Table  2.1  Threshold  energies  of  the  various  levels 

22 

of  Na  . 

At  three  different  bombarding  energies  <,  the  neutrons 

22 

feeding  the  various  levels  in  Na  have  the  following  energies 
at  0°  . 


Bombarding 

Energy 

Neutron  energy 

at  0°  for  the  various  levels 

ground  state 

587  Kev 

892  kev 

4.5  MeV 

2.20  MeV 

1.557  MeV 

1.211  MeV 

5  MeV 

2.687  MeV 

2.044  MeV 

1.704  MeV 

5.5  MeV 

Not  computed 

_ l _ 

2.192  MeV 

Bombarding 

Energy 

Neutron  energy  at  0°  for  the  various  levels 

1.532  MeV 

1 . 94  MeV 

1.95  MeV 

I.988  MeV 

4.5  MeV 

5  MeV 

5.5  MeV 

0.436  MeV 

0.966  MeV 

1.470  MeV 

No 

0.453  MeV 

0.989  MeV 

t  excite 

0.442  MeV 

0.979  MeV 

d 

0.390  MeV 

0.933  MeV 

Table  2.2  Energy  of  the  neutrons  feeding  the  various  levels 

22  o 

in  Na  and  emitted  at  0  from  the  beam  axis. 
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c)  Neutron  gamma  discrimination 

The  neutron  pulse  shape  discriminator  was  first 

tested  with  plutonium  -  beryllium  and  Co^  sources.  With 

the  discriminator  level  set  at  0,  it  was  checked  that  the 

rejection  ratio  was  as  high  as  claimed  by  the  Nuclear 

Enterprises  specifications  (500:1).  The  results  of  this 

test  are  shown  in  fig.  2.6.  The  ungated  spectrum  corresponds 

to  the  linear  output  only,  the  gated  spectrum  to  the  linear 

output  in  coincidence  with  the  pulse  shape  discriminator 

output.  A  second  test  was  run  with  a  thorium  C  source  whose 

maximum  gamma  ray  energy  is  2.62  MeV.  It  showed  that  the 

discriminator  level  had  to  be  set  to  12  in  order  to  get  a 

comparable  rejection  ratio.  Since  the  most  energetic 

strong  gamma  ray  from  the  reaction  (F^  +  a  )  is  the  2.1  MeV 

gamma  ray  corresponding  to  a  transition  from  the  second 

22 

excited  state  in  Ne  ,  a  discriminator  level  value  of  9 
was  found  to  be  optimum  during  testing  with  the  beam.  The 
results  of  this  last  test  are  illustrated  by  the  fig.  2.7 
and  2.8. 

d)  Yield  curves  for  the  0.892  and  1.532  MeV  states. 

Yield  curves  were  plotted  for  the  0.892  and  1.532 

22 

MeV  gamma-rays  from  Na  from  4.5  to  5-4  MeV.  The  energies 
were  determined  from  nuclear  magnetic  resonance  frequencies. 


COUNTS  PER  CHANNEL 
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60 


Fig  2.6:  Gated  and  ungated  neutron  spectra  for  PuBe  and  Co 


Current:  0.5 
Cha rge  :  250  gcou  1  . 
Target  :  310  pg/cm 
Discriminator  :  9. 
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Ungated  neutron  spectrum 
Bombarding  energy:  4.72MeV. 
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The  corrected  energies  show  a  discrepancy  of  about  150  kev 
compared  to  the  generating  voltmeter  readings.  Since  no 
other  monitoring  but  the  integrated  charge  on  the  target  is 
possible,  runs  were  taken  from  4.5  to  5.4  MeV  increasing 
by  steps  of  approximately  20  kev  and  for  an  integrated 
charge  of  50  microcoulombs.  After  these  runB  were  completed 
a  few  runs  were  again  taken  for  intermediate  energies  to 
check  the  reproducibility  of  the  data. 

o 

The  CaF2  loading  was  about  100  ^ig/cm  thick 
deposited  on  a  0.002  inch  thick  platinum  backing  and  the 
gamma-radiation  was  detected  at  90°  to  the  beam  direction. 

It  is  interesting  to  compare  these  yield  curves  to  those 
obtained  by  Freeman  and  Mani  (op.  cit),  in  the  region  4.5  to 
4.9  MeV.  Sinoe  the  target  used  by  Freeman  and  Mani  was 
about  five  times  thinner  than  the  target  used  in  this 
experiment,  one  should  expect  a  broadening  of  the  resonances 
in  the  present  investigation,  It  can  be  seen  that  the  two 
graphs  show  some  agreement  for  the  S92  kev  state  (fig,  2,9) 
but  do  not  compare  for  the  1,53®  MeV  state  (fig,  2,10), 

It  must  be  pointed  out  that,  in  the  present  experiment  an 
approximation  was  made  in  the  calculations  1  a  linear  back¬ 
ground  was  assumed  and  subtracted  from  the  total  number  of 
counts  in  the  peak,  Since  the  1,532  MeV  level  is  not 

i 

strongly  excited,  this  method  may  have  partially  affected 


o  • 
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Fig.  2.10:  Yield  of  the  1.532  MeV  transition. 
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the  results  whereas,  in  the  experiment  performed  by  Freeman 
and  Mani ,  the  computations  were  done  by  a  computer  with  a 
Gaussian  fit  program,  which  should  lead  to  more  reliable 
results  for  weakly  excited  levels. 

It  was  finally  decided  to  investigate  the  587 
and  892  kev  states  at  a  bombarding  energy  of  5*34  MeV 
corresponding  to  the  maximum  excitation  of  the  892  kev  level 
in  the  region  which  had  been  studied. 

e)  Angular  distributions. 

For  a  given  angle  a  singles  (non-coincidence) 
spectrum  was  taken,  then  a  coincidence  spectrum.  The  gamma- 
detector  whose  target-crystal  distance  waB  fixed  at  20.00  cm, 
was  then  moved  to  90°  and  a  singles  spectrum  was  taken  (as 
explained  in  section  2. 3. a).  The  target  chamber-detector 
distance  and  height  of  the  detector  were  checked  after 
each  position  change,  The  integrated  charge  was  fixed  at 
230  micro coulombs  and  runs  lasted  approximately  2  minutes, 

The  neutron  counter,  fixed  at  0°,  served  as  a  monitor  and  as 
the  coincident  neutron  detector,  Typical  singles  spectra 
for  a  bombarding  energy  of  5.38  MeV  are  shown  in  fig.  2,11 
and  2.12,  At  this  energy  the  triplet  situated  at  approxi¬ 
mately  1.95  MeV  is  excited  but  the  tables  of  Endt  and 
van  der  Leun  (op.  cit.)  show  that  gamma-rays  from  this  state 
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Fig  2.11:  Ungated  spectrum  for  a  bombarding  energy  of  5.3®  MeV. 
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cascade  to  the  ground  via  the  660  kev  state.  The  transition 

from  the  triplet  to  the  660  kev  state  corresponds  to  an 

energy  of  about  1 . 290  MeV.  This  transition  cannot  be  seen 

because  of  the  strong  competing  transition  of  the  1.28  MeV 

22 

level  from  Ne  .  The  illustration  2.12  corresponds  to  the 
same  spectrum  but  with  a  lower  gain  setting.  It  shows  the 
two  high  energy  gamma  rays  from  the  (a,p)  reaction  at 
2.1  MeV  and  3-35  MeV. 


f)  Angular  correlations  of  the  587  and  892  kev  states. 

The  axis  of  the  neutron  counter  was  at  90°  from 
the  beam  axis  and  at  a  distance  of  7»70  cm  from  the  target 
axis.  This  destroys  axial  symmetry  but  allows  runs  at  small 
gamma  detector  angles  to  be  performed,,  and,  as  seen  earlier, 
the  fundamental  assumptions  of  the  angular  correlation  theory 
can  still  be  considered  as  valid.  A  typical  spectrum  as 
shown  in  fig.  2.13  would  correspond  to  an  integrated  charge 
of  250  microcoulombs  and  a  running  time  of  approximately 
20  minutes. 

Remarks : 

l)  The  illustration  2.14  shows  a  background  spectrum 
taken  after  the  target  had  been  used  for  about  24  hours. 

Three  peaks  can  be  seen  which  correspond  to  potassium  40  (1.46  MeV), 
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Fig  2.13:  Neutron  gated  spectrum  for  a  bombarding  energy  of  5.38MeV 
3"x3"  gamma-detector  at  90  from  the  beam  axis. 
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sodium  22  (1.28  MeV)  and  (511  kev) .  The  sodium  peaks  had 
been  induced  from  the  reaction  as  was  checked  by  placing  a 
piece  of  lead  1  inch  thick  between  the  target  chamber  and 
counter.  The  potassium  40  activity  comes  from  the  environ¬ 
ment  and  the  detector  itself. 

2)  Two  points  have  to  be  mentioned  fion  the  study  of 
the  spectra: 

i)  From  the  figure  2.12,  two  high  energy  gamma-rays 

can  be  seen:  the  first  one,  at  3*35  MeV,  is  assumed  to  be 

22 

the  transition  from  the  second  excited  state  of  Ne  to  the 

ground  state.  This  transition  is  shown  in  Endt  and  Van  der  Leun 

tables  (En  62)  but  does  not  appear  in  more  recent  work  by 

Eswaran  and  Broude  (Es  64).  The  second  gamma-ray,  which  is 

assumed  to  be  due  to  a  transition  from  the  2.2  MeV  state  of 
22 

Na  ,  is  not  mentioned  in  the  literature.  However,  there  is 
not  enough  experimental  evidence  to  verify  these  assumptions. 

ii)  In  the  gated  spectra  a  peak  is  present  which 

corresponds  to  an  energy  of  220  kev  (fig.  2.4  and  2.5). 

This  transition  might  correspond  to  a  decay  from  the  2.2 

22 

MeV  level  in  Na  to  one  of  the  levels  of  the  triplet.  In 
this  case  again,  not  enough  experimental  evidence  is  avail¬ 
able  to  support  this  assumption.  Moreover,  in  the  second 
series  of  runs,  the  yield  was  too  low  to  make  a  further 
investigation  possible. 
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2.4  Calculations  and  Corrections. 

a.  Area  of  the  peaks. 

The  areas  of  the  peaks  were  computed  by  the  S.D.S. 

920  computer  from  a  Gaussian  fit  program  written  by  W.G.  Davies 
at  the  Nuclear  Research  Centre.  The  assumed  nature  of  the 
function  which  fits  the  background  is  entered  with  the  input 
data  and  may  be  linear.,  quadratic,  exponential  or  a  combination 
of  these.  After  a  few  trials  it  appeared  that  best  fits  were 
obtained  with  a  quadratic  function  describing  the  background. 
The  normalization  with  respect  to  the  neutron  monitor  was 
then  applied.  This  normalization  is  tagged  as  "standard 
normalization."  Independently,  a  second  correction  was 
applied  which  will  be  referred  to  as  "correction  for  variations 
at  90°"  (see  chapter  III). 

b.  Correction  for  neutron  counter. 

As  seen  earlier,  there  is  some  absorption  caused 
by  the  neutron  counter  for  small  angle  positions.  The  absorp¬ 
tion  coefficients  were  determined  by  measuring  the  ratio  of 
the  intensities  of  the  587  kev  and  892  kev  peaks  for  neutron 
counter  in  and  out,  at  a  given  position.  The  validity  of  the 
results  shown  in  the  next  table,  was  verified  by  doing  the 
same  measurements  with  Na22  and  Co^°  sources  and  checking 
the  continuity  of  the  data. 
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Absorption 

coef fi cient 

at 

0° 

15° 

0 

0 

on 

;oa-  for  i) 

587  kev 

1.20 

1 — 1 

0 

1 — 1 

1.00 

ii) 

892  kev 

LA 

1 — 1 

0 

1 — 1 

00 

0 

0 

1 — 1 

1.00 

Table  2.3  Correction  factors  for  absorption  in  the 

neutron  counter. 


c.  Correction  for  target  thickness. 

During  the  experiments  the  target  was  either  at 
45°  or  135°  with  respect  to  the  beam  axis. 

For  each  positions  the  gamma-rays  which  are  detected 
traverse  a  different  thickness  of  the  backing  and  this  must 
be  taken  into  account  when  corrections  are  applied.  Details 
of  calculations  will  be  given  in  the  appendix  A  2  and  results 
as  follows  for  the  two  investigated  gamma  rays. 


Angle 

of  detection  defined  with 

respect 

to  the 

normal 

to  the  target  (a) 

0 

y-ray 

energy 

0° 

15° 

30° 

45° 

O 

O 

VO 

75° 

587  kev 

0.997 

0.997 

0.9965 

0.996 

0.993 

0.982 

892  kev 

0.998 

0.998 

0.998 

0.9975 

0.995 

0.985 

Table  2.4.  Target  thickness  corrections  for  a  0.002” 

thick  platinum  backing.  Gamma  rays  are 
detected  at  20  cm  from  target. 
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d.  Angular  correlation  and  angular  distribution  fits. 

The  computation  of  the  Legendre  polynomial  co¬ 
efficients,,  A^.,  was  carried  out  by  the  nuclear  research 
centre’s  LGP-30  computer  from  a  program  written  by 
W.  G.  Davies.  This  program  will  accept  any  positive  angle 
from  0°  to  l80°  and  will  fit  the  data  to  any  set  of  five 
A^Ts  from  Aq  to  A^  providing  that  there  is  one  more  data 
point  than  A^’s.  For  the  analysis  of  experimental  data, 
the  fit  was  limited  to  fourth  degree  even  Legendre  poly¬ 
nomial  coefficients,  i.e.  A  ,  A^  and  A^ .  For  each  set  of 
data,  two  different  kinds  of  analysis  were  performed: 

i)  analysis  of  the  data  corresponding  to  the  ten 
measured  positions  from  0°  to  135°  by  steps  of  15°. 

ii)  analysis  of  the  data  corresponding  to  a  set  of 
seven  angles  from  0°  to  90°  and  in  which  the  data  corres¬ 
ponding  to  positions  symmetrical  with  respect  to  90°  had 
been  averaged, 

The  output  data  are  as  follows: 

p 

i)  best  X  fit  for  the  experimental  points 
ii)  unormalized  coefficients  AQ ,  Ag  and  A^ 

.  ?  t 

iii)  coefficients  Ag  and  A^  normalized  with  respect  to  AQ. 

All  the  coefficients  are  given  with  their  corresponding 

statistical  and  total  errors. 
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e)  Solid  angle  correction  for  gamma-ray  detectors. 
In  section  1.2.b,  it  has  been  seen  that  the 
angular  distribution  may  be  expressed  by: 


W (6)  =  1  +  a  2^2 ( cos  0)  +  a  ;,P;,(cos  0)  + 


4  4 


where  a  2  and  a  4  are  normalized  with  respect  to  a 


o 


The  a^  s  coefficients  are  related  to  the  expert 


mental  coefficients,  A  ^s,  by 


ak  \  =  Ak 


where  is  the  attenuation  factor  for  the  counter  which 

detects  the  gamma  rays  and  is  defined  by: 


J 

Qk  =  J 


k 


and 


rW 

Jk=  G  (£)  Pk  (cos  i)  sin  i  d  £ 

J  o 


e(£)  is  the  efficiency  of  the  counter  for  a  gamma-ray 
propagating  in  a  direction  at  angle  £  to  the  axis  of  the 


counter. 
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Davisson  (Da  64)  has  calculated  and  tabulated 
these  coef fi cients ,  for  a  3"  x  3"  detector.  But  since 

these  tables  are  limited  to  a  maximum  source  to  crystal  face 

T 

distance  of  15  cm,  the  values  of  s  had  to  be  extrapolated 
for  a  distance  of  20  cm. 

Values  of  are  equal  for  587  kev  and  892  kev 
gamma-rays  and  are  given  by: 

Qg  =  0.987 
=  0 . 950 
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CHAPTER  III. 


Experimental  results. 

a.  Validity  of  the  different  ways  of  analyzing  the  data. 

Following  the  methods  described  in  Chapter  II, 
the  angular  distributions  and  the  angular  correlations  have 
been  analyzed.  Figures  3*3  to  3*10  (inclusive)  show  the 
different  angular  correlation  and  distribution  data  and  the 
calculated  best  fits. 

It  is  interesting  to  compare  the  different  ways 
the  data  have  been  analyzed  for  a  given  transition,  namely: 
i)  by  considering  the  neutron  monitor  normalization 

("standard  normalization")  first  fitting  on  10  points 
and  later  on  7  points  only. 

ii)  by  applying  a  second  correction  monitored  by  the 
relative  intensity  of  the  peak  from  the  ungated 
spectrum  taken  at  90°  from  the  beam  axis  ("correction 
for  variations  at  90°"). 

The  following  two  illustrations  (fig.  3*1  and  3*2) 
give  the  results  of  the  computations.  The  "a"  coefficients 
are  plotted  with  their  errors  corresponding  to  the  larger 
of  statistical  and  total  errors  . 

*Note :  the  total  error  as  defined  by  Rose  (Ro  53)  is  a  com¬ 
parison  between  the  true  coefficient  a^  and  the  least  squares 
value  a^.  This  error  is  shown  to  be  defined  by  the  product 
aX  where  cr(statistlcal  error)  is  the  root  mean  square 
deviation  of  the  a’  . 

D 


•'  tl  ‘ 


angular  distribution  angular  correlation 
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Comparison  of  anqular  distribution  and  angular 
correlation  coefficients  for  the  587  keV  state. 
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Comparison  of  angular  distribution  and  angular 
correlation  coefficients  for  the  892  keV  state, 
(cf.  previous  plot) 
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An  important  parameter  to  consider  first  is  the 

2 

best  fit  coefficient,  "X  .  This  coefficient  describes  the 
closeness  of  the  experimental  curve  to  a  function  of  the 
fourth  order  of  even  Legendre  polynomials  (i.e., 

W(0)  =  1  +  a  2?2  (cos  0)  +  a  cos  0)  and  is  given  by 

the  minimum  of  the  expression: 

4 

2  W1(M,1  -  2  a^P^  (cos  0i)  )2 

i  ^even 

t  i 

where  the  |_i.  s  are  the  observed  values  and  w.  s  the 
i  i 

corresponding  weights. 

2 

As  pointed  out  by  Rose  (Ro  53) *  a  value  of  X 

appreciably  greater  than  one  implies  either  systematic 

errors  or  the  use  of  an  incorrect  function  in  the  least 

2 

squares  fit.  Conversely,  a  value  of  *X  much  smaller  than 
one  indicates  either  that  too  many  polynomials  are  intro¬ 
duced  to  fit  the  experimental  graph  or  that  one  of  the 
experimental  data  has  an  error  much  smaller  than  the  others. 
(This  smallest  error  is  thus  taken  as  the  most  significant 
weight ) . 


As  seen  from  fig.  3*1  and  3.2,  the  only  cases  of 

poor  fits  are  those  of  the  angular  distributions.  Except 

2 

for  one  case,  the  X  values  are  all  smaller  than  one  for 
the  angular  correlations.  Because  statistics  are  better 
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for  the  distribution  than  the  correlation  measurements ,  some 
systematic  errors,  such  as  assymetry  of  the  correlation  table, 
are  much  more  revealed  by  distribution  data  and  are  the 
cause  of  this  difference  in  the  best  fit  coefficients. 

Secondly,  it  is  interesting  to  see  how  much 
improvement  has  been  obtained  by  the  second  correction. 

From  the  various  plots  of  angular  distributions  and  from 
fig.  3.1  and  3*2,  one  can  see  that  the  second  correction  led 
to  a  better  fit  for  the  two  angular  distributions  and  for 
the  587  kev  angular  correlation,  although  the  improvement 
in  the  latter  case  is  not  as  large.  On  the  other  hand,  the 
892  kev  angular  correlation  fit  got  worse  after  the  second 
correction  had  been  applied.  This  second  correction  has 
been  applied  because  it  was  noticed  that  the  relative 
intensities  of  the  587  and  892  kev  peaks  from  non-coincidence 
spectra  taken  at  a  position  of  90°  from  the  beam  axis,  were 
not  constant.  The  complete  correlation  experiment  was 
performed  in  three  days.  For  the  longest  day's  run,  the 
maximum  absolute  variation  of  the  892  kev  peak  relative 
intensity  reached  13$.  However,  it  must  be  pointed  out 
that  the  relative  intensity  of  this  peak  was  defined  with  an 
average  relative  error  of  4$.  This  correction  is  used  as 
a  normalization  method  in  angular  correlation  experiments 
dealing  with  natural  radioactive  sources  because  finite 
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lifetimes  are  then  involved.  For  induced  reactions,  where 
machine  time  may  be  limited,  this  method  could  hardly  be 
applied  if  long  runs  are  required  for  good  statistics. 
Moreover,  since  the  angle  positions  are  taken  at  random  this 
kind  of  error  should  cancel  out  after  averaging  all  the  data 
corresponding  to  one  position. 

As  a  conclusion  for  this  first  part,  it  can  be 
said  that  the  normalization  with  respect  to  the  neutron 
monitor  and  the  usual  corrections  (target  thickness,  finite 
solid  angle  of  detector)  are  necessary  and  sufficient. 

A  second  point  which  would  have  proved  interesting 
to  investigate  was  whether  the  results  from  the  measurements 
of  the  587  kev  transition  lead  to  the  known  1+  assignment. 

Unfortunately,  such  a  conclusion  cannot  be  reached 
because  of  three  reasons: 

i)  large  errors  of  the  experimental  results 

ii)  multiple  feeding  of  the  587  kev  level 

iii)  long  lifetime  of  the  66 0  kev  and  587  kev  levels. 

These  three  causes  will  now  be  examined. 

i)  Large  errors  of  the  experimental  results. 

Fig  3.3,  and  3»5  show  two  different  computed 

fits  for  the  experimental  data:  one  is  limited  to  the  fourth 
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order  of  even  Legendre  polynomials ,  the  latter  to  the  second 

order  only.  Because  of  the  large  errors,  the  value  of  the 
2 

X  best  fit  coefficient  remained  either  much  larger  than  one 
(for  angular  distribution)  or  much  smaller  than  one  (for 
angular  correlation).  From  these  different  computations, 
the  values  of  the  a’  coefficients  are  given  with  relative 
errors  ranging  from  20  to  120$.  Such  large  errors  would 
probably  lead  to  ambiguous  spin  assignment  even  if  no  other 
corrections  were  to  be  applied. 

ii)  Multiple  feeding  of  the  587  kev  level. 

It  can  be  seen  from  diagram  2.1  that  the  587  kev 
level  is  not  only  fed  by  neutrons  but  also  by  gamma-rays 
cascading  via  the  660  kev  level  from  the  2.2  MeV  state  and 
the  triplet  at  1.9  MeV.  One  should  then  consider  the 
complete  correlation  and  distribution  formulae  in  which  the 
unobserved  transitions  are  taken  into  account.  As  shown  by 
Devons  and  Goldfarb  (De  57) >  the  formulae  can  still  be 
reduced  to  the  form  W(0)  =  2  a^Q^P^cos  0)  but  then,  the 
a^.s  also  describe  the  unobserved  radiations  through  various 
Racah  and  Clebsh-Gordon  coefficients.  It  has  been  seen 
earlier  (Chapter  i)  that  if  a  state  fed  only  by  neutrons  has 
spin  0  or  1,  the  coefficients  of  angular  distributions  and 
angular  correlations  are  identical  because  no  further  restriction 


b  • 
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on  the  number  of  magnetic  substates  which  are  fed  is  brought 
by  the  correlation.  Since  the  2.2  MeV  level  and  the  triplet 
levels  may  have  spins  different  from  0  and  1,  one  should 
then  expect  different  a  coefficients  for  angular  distribu¬ 
tion  and  angular  correlation.  However,  for  a  bombarding 
energy  of  5.38  MeV,  these  levels  are  probably  weakly  excited, 
(it  must  be  recalled  that,  in  the  last  series  of  runs,  no 
transition  corresponding  to  an  energy  of  more  than  892  kev 
was  observed  on  the  n-y  coincidence  spectra  because  the 
yield  was  too  low  for  the  running  time).  Consequently, 
some  assumptions  might  be  made  if  this  were  the  only  correc¬ 
tion  to  apply.  However,  as  will  now  be  seen  another  per¬ 
turbing  condition  is  happening. 

iii)  Long  lifetimes. 

According  to  Temmer  and  Heydenburg  (Te  58),  the 

660  kev  level  has  a  half-life  Ti  equal  to  14  nanoseconds 

2 

and  the  587  kev  a  Ti  equal  to  266  nanoseconds.  As  seen 

2 

earlier  (section  2.1. b)  the  resolving  time  of  the  coincidence 
circuit  is  approximately  40  nanoseconds.  As  a  result  of  the 
lifetime  of  the  587  kev  excited  state  being  much  longer 
than  the  coincidence  resolving  time,  the  n-y  coincidence 
counting  rate  will  be  lowered.  As  a  matter  of  fact,  the 
problem  is  more  complex  because,  as  shown  by  Steffen  and 
Frauenfelder  (St  64),  perturbations  of  the  correlations  are 
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possible  when  the  lifetime  of  the  intermediate  state  is  longer 
than  10  seconds.  Perturbed  correlations  have  been  studied 
extensively  in  the  recent  years  and  the  most  relevant  papers 
on  this  topic  can  be  found  in  (Ka  64).  However*  these 
studies  have  been  mainly  directed  towards  the  study  of 
radioisotopes  and  the  influence  of  their  physical  and 
chemical  environment. 

If  one  considers  a  cascade  I.  -►  I  -*•  If*  the 

angular  correlation  of  gamma  rays  will  be  altered  if  the 

nuclei  in  their  intermediate  state  are  subject  to  torque* 

due  to  the  interaction  of  either  the  magnetic  dipole  moment 

with  an  extranuclear  field  B  or  the  electric  quadrupole 

62V 

moment  Q  with  an  electric  field  gradient  — 75-  .  These 

6Z2 

interactions  produce  a  precession  of  the  nuclei  around  the 

symmetry  axis  which  results  in  a  time  dependence  in  the 

density  matrix  of  the  intermediate  state.  In  solids*  the 

magnetic  interaction  can  usually  be  neglected  in  comparison 

with  the  quadrupole  interaction  (St  64).  The  magnitude  of 

the  quadrupole  interaction  is  described  by  the  precession 

62V 

frequency  ^  which  is  proportional  to  Q  and  — «• 

6Z 

The  correlation  function  is  multiplied  by  an 
attenuation  factor  Gk(f)  which  is  a  function  of  the  parameters 
of  the  states  involved  and  the  precession  frequency.  Generally* 
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this  attenuation  factor  is  determined  experimentally 
(interesting  examples  are  given  in  (Fr  55)  )  and  leads  to 
the  value  of  the  electric  quadrupole  moment  of  the  state 
which  is  studied. 

In  the  present  case,  the  situation  is  even  more 
complex  since  two  long  lifetime  states  are  involved 
and  although  it  might  be  possible  to  assume  some  value  for 
the  electric  field  gradient  due  to  surrounding  electrons, 
the  electric  quadrupole  moments  are  not  known. 

In  conclusion,  very  little  can  be  said  about  a 
possible  assignment  of  the  587  kev  level  from  the  experi¬ 
mental  results.  Long  calculations  and  many  simplifying 
assumptions  would  be  necessary  and  are  not  worthwhile 
because  of  the  large  errors  of  these  results. 

b.  Spin  assignment  of  the  892  kev  level. 

The  situation  for  this  level  is  fortunately  clear 
cut:  both  values  of  a'2  and  a’^  are  well  defined,  no 

multiple  feeding  is  involved.  A  question  may  be  raised 
however  concerning  the  half  life  of  this  level.  No  datum 
is  actually  given  but  compared  with  other  nuclei  [a  table  in 
(Ka  64)  gives  all  the  measured  lifetimes  for  the  levels  of 
all  nuclei]  it  is  a  good  assumption  that  the  half  life  of  the 
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corrections  for  variations  at  90  degrees 
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892 


kev  level  is  less  than  10 


-11 


sec. 


From  the  best  fit  corresponding  to  the  analysis 
of  7  angles  using  the  "standard  normalization",  the  values 
of  the  Legendre  polynomial  coefficients  are  given  by: 


a ' 2  =  0.4l8  ±  0.047 
a 1 =  0.223  ±  0.048 


('X.2  »  0.69) 


This  angular  correlation  is  represented  in  graph 
3.8.  The  continuous  line  represents  the  computed  best  fit. 

The  dotted  line  corresponds  to  a  best  fit  limited  to  the 
second  degree  Legendre  polynomial  P^fcos  0).  This  second 

4- 

fit  is  very  poor  and  thus  eliminates  a  priori  a  1  assignment. 

? 

From  the  sign  of  a  two  different  assignments 
are  possible:  2  or  4  ,  Qualitatively,  the  2  is  more 
probable  because  there  is  a  small  discrepancy  between  the 
Legendre  polynomial  coefficients  for  both  angular  distribu¬ 
tion  and  angular  correlation  (cf.  fig.  3*2).  Better  informa- 

f  t 

tion  can  be  obtained  from  the  plots  of  a  ^  versus  a  ^  where 
2 

X  is  fixed  and  5  is  taken  as  the  variable.  For  these 

2 

plots,  a  few  values  of  the  population  parameter  ratio,  X  , 
have  been  considered  including  the  two  extreme  values 
X2  =  0  and  X2  =  °°  ,  (fig.  3.11  and  3.12).  From  graph  3*11 
(transition  4  to  3)j  it  can  be  seen  that  the  area  representing 


ig  3.11:  Plots  of  a'2  versus  aA4  for  all  values  of  6 
and  for  a  4  to  3  spin  change,  assuming 
s-wave  neutrons  only. 

Different  values  of  the  population  para¬ 
meters  ratio,  x  ,are  considered. 
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Fiq  3  12:  Plots  of  e'2  versus  a^for  all  values  of  6  and  for 
a  2  to  3  spin  change  assuming  s-wave  neutrons  only. 

Different  values  of  the  population  parameters  are 
cons i de  red . 
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0.1  1  10 


Fig  3.13:  A  plot  of  x2  versus  o  using  the  measured  values  of 
a!,  and  aj  for  the  892keV  gamma  ray  and  the  predic 
ted  a  coefficients  for  a  2  to  3  spin  change. 
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t  » 

the  experimental  values  of  a  ^  and  a  4  does  not  correspond 

to  any  physical  solution  since  it  does  not  intersect  the 

curves.  From  graph  3*12  (transition  2  to  3)*  one  can  see 
2 

that  for  X  =  0  there  is  an  intersection  corresponding  to 
an  approximate  value  of  6  equal  to  0.5*  One  can  then  con¬ 
clude  that  the  spin  2+  is  the  only  possible  assignment  for 
this  state. 

Because  there  are  two  unquantized  parameters 

involved  in  the  correlation^  the  mixing  ratio  5  is  not 

2 

sharply  defined  The  figure  3° 13  gives  the  plot  of  X 

T  T 

versus  5  for  the  two  experimental  values  of  a  2  and  a  ^ . 

The  measured  values  a  2  and  a  ^  are  put  in  the  Legendre 

polynomial  equations  for  the  transition  2  to  3  (see  Appendix 

A  l)  and  X  is  then  determined  for  a  proper  range  of  5. 

Two  regions  are  thus  defined  and  their  intersection  gives 

2 

the  most  probable  range  of  X  and  5.  These  values  are  given 
by: 

x2  =  0  to  0.40 

5  =  0. 4l  to  1.03  or  0.72  ±  0.31 


They  will  now  be  discussed. 
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c.  Discussion. 

By  considering  the  possible  value  of  the  ratio  of 

the  population  parameters  P(±l)  and  P(0),  it  may  be  assumed 
2 

that  X  =  0  is  not  likely,  because  this  would  correspond  to 

the  0  substate  being  populated  only.  It  must  be  recalled 

that  all  the  calculations  suppose  that  the  neutron  counter 

is  at  0°  and  that  only  S-wave  neutrons  are  considered,, 

thus  limiting  the  magnetic  substates  populated  to  0  and  ±1. 

Actually,  because  of  the  finite  size  of  the  detector,  the 

±2  magnetic  substates  may  be  weakly  populated.  However, 

calculations  show  that  the  probability  of  populating  these 

substates  is  small,  a  few  percents  generally.  Assuming  that 

2  2 

the  mean  value  of  X  ,  X  =  0.2,  is  more  probable  than  the 
others  and  supposing  that  the  ±2  magnetic  substates  are  not 
populated,  one  can  calculate  the  population  parameters  of 
the  0  and  ±1  substates.  Since  2  P(l)  +  P(0)  =  1,  one  gets 
P(0)  =  71$  and  P(l)  =  14.5$.,  if  X2  =  =  0.2.  This 

result  is  reasonable  compared  to  the  equivalent  data 
obtained  from  other  experiments  (see  e.g.  (Li  6o)  ). 

The  value  of  8  which  is  the  quadrupole  to  dipole 
mixing  ratio  is  given  by  0.72  ±  31  which  corresponds  to  a 
relative  error  of  43$.  As  stated  earlier  (section  I.3,c), 
Morpurgo  (Mo  63)  has  considered  the  Weisskopf  estimates  and 
showed  that,  in  the  case  of  self -mirrored  nuclei,  the 


87 


magnetic  dipole  transition  with  AT  =  0  should  be  reduced  by 

a  factor  150.  It  is  also  well  known  that  compared  to  the 

Weisskopf  estimates  the  electric  quadrupole  is  generally 

enhanced.  For  the  892  kev  transition,  the  Weisskopf  esti- 

.  _4 

mate  predicts  a  ratio  E2/M1  3  x  10  .  Despite  its  in¬ 

accuracy,  the  present  result  agrees  with  the  fact  stressed 
by  Morpurgo  that  one  should  expect  a  strong  Ml  -  E2  mixing 
even  for  rather  low  energy  transitions. 


The  2+  spin  assignment  suggests  that  the  892  kev 

state  is  the  first  level  of  the  K  =  2+  rotational  band  or 

+  + 

the  second  level  of  the  K  =  1  band.  The  K  =  2  rotational 
band  is  built  on  the  transition  of  one  odd-nucleon  from 
the  3/2  +  Nilsson  level  to  the  1/2  +  level.  Using  a  value 
of  Tjfj  =  0.23  MeV,  as  proposed  by  Paul  (op.cit.)  would  rule 

out  the  assumption  that  the  892  kev  level  is  the  second  level 

+  fi2 

of  the  K  =  1  rotational  band.  This  value  of  is  similar 


20  24 

to  the  value  for  Ne  and  Mg 


2  J 
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Ar  65 
Br  60 
Da  52 
Da  64 
De  57 
En  62 
Es  64 
Fr  64 
Go  60 
Gov  60 
Ka  64 
Li  60 
Li  6l 
Li  64 
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APPENDIX  A-l 


Angular  correlation  formulae  assuming  s-wave  neutrons, 

for  different  possible  spin  assignments,  in  the  case  of  the 

1 Q  22 

reaction  P  ^(any)Na 


5  is  the  mixing  ratio. 

2 

x  the  ratio  of  the  population  parameter  of  the 
magnetic  sublevels  ±1  over  the  population  parameter  of  the 
sublevel  0:  |x2  =  . 

It  is  reminded  that  the  ground  state  is  a  3+  and 
that  the  parity  of  the  unknown  states  involved  is  positive. 


Spin  1+  : 


w(e) 


,  ,  0.144  +  1.170  5  ■ 

1  +  . .  . . .  . K' 

1+5^ 


0.250  6*  l-x‘ 


l+2x 


2  P2(cos  0) 


Spin  2+  s 


W(6)  «  1  - 

1+6" 


l+x‘ 


l+2x 


2  pg(oos  e) 


+ 


1.219  -  1.626  xs 

(l+2xa)  "(1+B2) 


(cos  B) 


Spin  3+  ‘ 


A  2 


W(0)  i  +  (°»  500-0 » 998  5  +  0.111  52X0. 748-1.496  5-0.168  s2)x: 

(1  +  2x2)  (1  +  62) 


P0(cos  0)  -  (0-53.3  +  0.180  x  ) 


1  +  2x 


1+6 


2  P/|  (cos  0) 


Spin  4+  : 

W(0)  =  1  + 


-(0.537+0. 6o8x2 ) - ( 2 . 138+3 • 632x2 ) 6+ ( 0 . 051+0 . 086x2 ) 6 2 

o  o 


(1  +  2x  )  (1  +  6  j 


Po  (cos  e)  +  0-732(l+x  ) 


75-  - p  Pk(cos  0) 
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The  following  five  graphs  represent  the  plots  of 
some  of  the  coefficients  of  the  Legendre  polynomials  versus 
the  mixing  ratio  for  different  population  parameters  ratios. 
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Fig.AJ.h  cu  /(L  vs 5  for  a  transition  I  to  3  ( s-wave  neutrons  only) 
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Fjg:  A. 1.2  :  o^y/c^  vs  S  for  a  transition  2  to  3 
(s-wave  neutrons  only) 
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Fig.  A. 1.3  :  £>2/^0  vs  ^  f°r  0  ♦rons»ti°n  2  to  3 
(s-wave  neutrons  only) 
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Fig.  A. 1.4  :  a4/6Q 


for  a  transition  3  to  3 
(swave  neutrons  only) 
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Fig.  A. 1.5: 


a  /  vs  <S  for  a  transition  4  to  3 
4/  0 

(s-wave  neutrons  only) 
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APPENDIX  A-2 

Corrections  for  target  thickness 


The  coefficients  of  correction  for  absorption  due 
to  target  thickness  were  calculated  for  a  platinum  backing 
0.002”  thick  (12.7  x  10“3  mm). 


The  beam  is  supposed  to  be  well 
defined  in  0.  The  absorption  of 
absolute  intensity  I  has  been 
calculated  with  respect  to  the  angle  a 
which  is  defined  with  respect  to  the 
normal  in  0. 


After  passage  through  a  thickness  X  the  reduced 
intensity  I  is  given  by: 

» 


where  \±  is  the  absorption  coefficient  in  cm  per  atom  if 

T 

the  thickness  X  is  expressed  in  atoms  per  square  centimeter. 

2 

Thickness  in  atoms/cm  : 


X ’  -  ^  x  P  x  N 
A 

.  -4  \ 

where  X  is  the  thickness  in  cm  (12.7  x  10  cm) 

o 

p  is  the  density  of  platinum  (21.37  g  /cmJ) 

OO 

N  is  Avogadro’s  Number  (6.02  x  10  J) 

A  is  the  atomic  number  of  platinum  (195*09). 


These  data  yield  a  value  of  X’  equal  to 
8.36  x  lO1^  atoms/cm^. 

The  absorption  coefficients ,  |i!s,  were  determined 
by  interpolation  from  the  data  published  by  Davisson  (Da  52) 
and  are  given  by : 

\i  =  36.8  barns/atom  for  E  =  587  kev 
[i  -  24.5  barns/atom  for  =  892  kev. 

Since  for  each  angle  a  the  corresponding  thickness 

» 

Xa  is  given  by 

x'  =  x' 
a  cos  a 

it  was  then  possible  to  plot  the  ratio  i/I  versus  a. 

For  a  3"  x  3"  detector  whose  crystal  face  is 
situated  at  a  distance  of  20.00  cm  from  the  target  axis, 
the  half  angle  limited  by  the  solid  angle  has  a  value  of  10.8 

Then  for  a  varying  from  0  to  90°  by  steps  of  15°* 
integration  was  performed  from  a  -  10.8°  to  a  +  10.8°  in 
order  to  determine  the  average  target  thickness  correction. 
The  results  of  these  computations  have  been  given  in 
table  2.4  (section  2,4, o). 


